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The  accuracy  of  measurements  or  predictions  of  the  frequency  response 
of  a single  degree  of  freedom  system  depends  on  the  frequency  resolution  of 
the  analysis  technique.  If  the  ratio  of  analysis  bandwidth  to  the  3dB  band- 
width of  the  resonance  is  too  great,  negative  bias  errors  occur  at  resonance 


in  estimates  of  both  power  spectral  levels  and  transfer  function  amplitudes. 
These  errors,  which  have  been  calculated  numerically,  are  presented  and 
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INTRODUCTION 


The  accurate  measurement  and  prediction  of  the  behaviour  of  a single 
degree  of  freedom  system  is  important  in  any  field  in  which  the  behaviour 
of  more  complex  structures  may  be  synthesised  by  superposing  the  behaviour 
of  many  single  degree  of  freedom  systems. 

In  the  field  of  structural  dynamics,  it  is  usual  to  measure  the 
experimental  frequency  response  of  a given  system  by  using  a dual  channel 
Fast  Fourier  Transform  (FFT)  analyser  with  cross  channel  computation  facili- 
ties. Such  an  analyser  carries  out  a Discrete  Fourier  Transform  on  sampled 
time  history  data  and  outputs  a frequency  response  which  is  defined  at  a 
finite  number  of  equally  spaced  discrete  frequencies.  The  spacing  between 
adjacent  frequency  estimates  may  be  defined  as  the  nominal  analysis  bandwidth. 
Theoretical  predictions  of  the  behaviour  of  the  same  system  may  typically  be 
made  using  the  Finite  Element  (FE)  method  to  generate  the  forced  response  of 
the  system.  In  these  circumstances,  the  theoretical  frequency  response  is 
normally  calculated  at  a series  of  chosen,  discrete  frequencies  whose  spacing 
may  once  again  be  defined  locally  as  the  analysis  bandwidth.  Values  of  the 
theoretical  frequency  response  at  other,  intermediate  frequencies  are  usually 
determined  by  linear  interpolation. 

In  either  analysis,  full  resolution  of  the  frequency  response  is  only 
possible  if  the  analysis  bandwidth  Af  is  small  enough  in  the  close 
neighbourhood  of  resonance  in  which  the  rate  of  change  of  both  amplitude  and 
phase  of  the  frequency  response  will  be  maximised.  For  theoretical 
predictions,  errors  are  clearly  introduced  whenever  linear  interpolation  is 
applied  to  a frequency  response  that  is  insufficiently  defined  as  a result  of 
the  frequency  estimates  being  too  widely  spaced.  For  experimental  measure- 
ments of  broad  band  data,  Forlifer  (l)  and  Bendat  et  al  (2)  have  shown  that 
tolerable  errors  are  introduced  into  power  spectrum  measurements  of  a single 
degree  of  freedom  (SDOF)  system  as  long  as  Af  satisfies  the  criterion, 

Afa  ^ 0.25 

Af 

s 

where  Afg  is  the  3dB  bandwidth  of  the  SDOF  system  resonance. 

However,  as  will  be  shown  in  Section  2,  the  errors  due  to  analysis 
bandwidth  are  much  more  significant  for  measurements  of  transfer  function 
amplitude  than  they  are  for  power  spectrum  measurements.  Fortunately,  the 
presence  of  these  increased  errors  will  be  indicated  in  experimental  measure- 
ments by  a corresponding  loss  of  coherence.  Expected  errors  in  transfer 
function  phase  are  also  shown  in  Section  2.  In  Section  3»  experimental 
measurements  taken  from  an  electrical  SDOF  system  are  presented  to  substantiate 
the  numerical  results.  The  implications  of  bandwidth  induced  errors  to 
structural  analysis  are  discussed  in  Section  4. 


. DERIVATION  OF  BIAS  ERRORS  IN  POWER  SPECTRUM  AND  TRANSFER  FUNCTION 
ESTIMATES  FOR  VARYING  ANALYSIS  BANDWIDTHg. 

Consider  a SDOF  system  which  is  described  by  the  second  order  differen- 
tial equation 

y + 25wNy  + wN2y  = x(t) 

where  wN  is  the  system  natural  frequency 

is  the  viscous  damping  factor,  stated  a6  a fraction  of  critical 
damping 

y is  the  displacement 

and  x(t)  describes  an  externally  applied  forcing  function,  per  unit  mass. 

Y(w),  the  Fourier  transform  of  the  displacement  is  related  to  X(w),  the 
Fourier  transform  of  the  forcing  function  by  way  of  the  frequency  response  or 
system  transfer  function  which  may  be  defined  as  H(w)  such  that 


H(w) 

= Y(w) 

= a(w)  + ib(w) 

x^TT 

, 2 2x 

where 

a(w) 

(w^  - w^) 

(w  2 - 
'WN 

2\2  ^ . <2  2 2 
w ) + 4^  w wN 

and 

b(w) 

s 

-2^wwn 

(w  2 - 
^WN 

2v2  , ,-2  2 2 

w ) + 4 % w wN 

The  magnitude  of  H(w)  is  defined  by 

I H(w)|2  = H(w)H*W)  = a2(w)  + b2(w) 
Also  at  the  discrete  frequency  w 
= a2(w)  + b2(w) 


where  G (w)  and  G (w)  are  the  power  spectra  of  x(t)  and  y(t)  respectively, 
xx  yy 

The  phase  of  H(w)  is  defined  by  0(w)  where 

0(w)  = tan”1  / b(w)  \ 

\ a(w) ) 

In  addition,  if  the  resonant  frequency  wfi,  is  defined  as  the  frequency  at 
which  the  forced  response  has  maximum  amplitude,  then 

»B2  . «„2(l-2f> 

and  Wg,  the  3dB  bandwidth  of  the  system  resonance  is  given  by 
wg  = 2£wn 


G (w) 

yy,  . 

G~  ~ ( w ) 
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w. 


and 


G (w) 

n. 


G (w) 

XX 


■^r 


P P 

a (w)  + d (w)  dw. 


M,w  = w 


12 


H(w12)  and  hj/w12^  &re  al8°  related  by 


H(W12) 


lim  Hu(w._) 

M 12 


In  the  neighbourhood  of  resonance  defined  by  full 

resolution  of  H(w)  requires  that  the  measured  estimate  H^Cw^)  should 

define  closely  the  system  behaviour  at  resonance.  In  this  context,  it 
is  convenient  to  define  AP(  Aw,  C , wR),  AH(  Aw,  6,  wR)  and  A0(  Aw,  £,  wR) 

respectively  as  errors  at  resonance  in  the  power  spectrum  amplitude, 
transfer  function  amplitude  and  transfer  function  phase.  Thus 


AP(Aw,e,  wR)  = 101og10  j 

G (w) 

• 

G (w,) 

yy  2 

~i 

t 

i 

G (w) 

XX 

k,w«,12 

G . (w  ’> 
Lxx  Rj 

All(  Aw,  £ , Wp) 

= 101°eio 

IW'2 

|H(wr)  | 2 

A0(  Aw,  fe,  wR) 

= 0(wR)  - 0M(w12) 

where 

Aw  = wg  - 

and 

6 = Wr  ~ W1 

(0*fc«l) 

W2  " W1 

p 

For  any  particular  measurement,  the  expected  coherence  "5S  ( Aw,  £,  wR) 
is  defined  by 

*2(Aw,  £,  wR)  = 10*10 

where  ^ = AH(  Aw,  £,  wR)  - AP(  Aw,  6,  wR) 

However,  in  practical  circumstances,  £ is  effectively  random  and 
as_a  result_it  is  more  convenient  to  utilise  the  mean  expected  errors 
AP(Aw),  AH(  Aw)  and  A0(  Aw)  and  their  respective  variances 

o'  2 ( Aw) , o-2h(  Aw)  and  ( Aw).  Ap(Aw)  and  o-2  (Aw)  are  defined 
P M V p 
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respectively  by 


AP(AW) 


o'  ^ ( Aw ) 

P 


AP(  Aw,  £,  wR) 


( AP(  Aw) 


d£ 


AP(  Aw,  £,  wR))2 


dC 


Values  of  Ap(Aw),  AH(Aw)  and  A0(Aw)  have  been  evaluated  on 
a PDP  11/40  minicomputer  by  numerical  integration  and  are  presented  in 
Figures  1,  2 and  3 in  which  the  mean  expected  errors  and  their  standard 
deviations  are  plotted  against  (Aw/p^w^.  To  carry  out  the  integration 

the  bandwidth  Aw  has  been  subdivided  into  200  bands  and  500  evenly 
distributed  random  values  of  & have  been  used.  Using  these  data,  the 
integration  is  convergent  and  accurate  to  four  significant  figures.  The 
same  numerical  integration  has  also  been  used  to  calculate  the  mean 
expected  coherence  ^(Aw)  and  values  of  its  standard  deviation.  These 
data  are  shown  in  Figure  4. 


It  is  apparent  from  Figures  1,  2 and  3 that  the  mean  errors  are 
insensitive  to  the  damping  factor  % , provided  that  1*0.1. 

2.2  Bias  errors  in  theoretical  predictions. 


The  mean  errors  which  affect  the  results  of  finite  bandwidth 
experimental  analysis  also  apply  to  theoretical  predictions  in  which  a 
broad  band  forcing  function  is  utilised.  However,  in  the  field  of 
structural  dynamics,  a discrete  frequency  forcing  function  is  normally 
used,  especially  for  finite  element  calculations.  In  these  circumstances, 
the  predicted  response  of  a SD0F  system  will  differ  from  the  true  resonant 
response  unless  the  discrete  frequency  of  the  input  is  exactly  equal  to 
the  resonant  frequency  of  the  system.  In  this  case,  it  is  simple  to 
define  the  errors  which  occur  in  estimates  of  the  resonant  values  of 
power  spectrum  amplitude,  transfer  function  amplitude  and  transfer 
function  phase  as  APT(w),  AH,j,(w)  and  A0T(w)  where  w is  the 

excitation  frequency.  If  a sinusoidal  force  of  frequency  w is  input 
to  a SDOF  system  of  natural  frequency  w^,  it  may  be  shown  that 


where 


AHt(w)  = APt(w) 

v-j,  i - A* 

iWN  N ' 


It  may  also  be  shown  that 


A0t(w) 


tan 


-1 


101og10ot 


(0.25«<*«1) 


frj w»„  \ *t>"  -1  (i 


I 


•I 


l 
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Values  of  Ah^Jw),  Ap^Cw)  and  A0T(w)  are  shown  plotted  against 

Cw  - wfjV^Wjj  in  Figure  5.  From  Figure  5,  it  may  be  determined  that  if 

the  Forlifer  criterion  (l)  is  adopted  in  these  circumstances,  then 
amplitude  errors  are  restricted  to  0.3dB  although  phase  errors  may  be  as 
great  as  l4°. 

3.  COMPARISON  OF  NUMERICAL  RESULTS  WITH  EXPERIMENTAL  DATA  FOR  TWO  SDOF  SYSTEMS. 

In  order  to  substantiate  the  numerical  results  presented  in  Section  2, 
experimental  measurements  have  been  carried  out  on  two  separate  electrical  SDOF 
systems.  Each  of  these  resonators  comprised  a resistance,  inductance  and 
capacitor  as  shown  in  Figure  6.  The  measured  system  transfer  function  H(w) 
is  defined  as 

H(w)  = v0UT(w) 


The  two  resonators,  which  had  natural  frequencies  of  500Hz  and  4kHz,  had 
respective  damping  factors  of  0.25  and  0.05. 

A Spectral  Dynamics  SD36O  dual  channel  FFT  processor  was  used  to  measure 
H(w)  for  varying  analysis  bandwidths.  The  nominal  analysis  bandwidth  was 
increased  by  a factor  of  1.8  by  the  SD36CS  standard  Kaiser-Bessel  data  smoothing 
window. 

Experimental  data  from  these  tests  are  compared  with  the  numerical  data  in 
Figures  7-10.  Figures  7 and  8 show  respectively  the  transfer  function  errors 
that  were  obtained  using  the  4kHz  and  500Hz  resonators,  whereas  Figure  9 shows 
the  power  spectrum  errors  which  were  derived  simultaneously.  In  all  cases,  the 
theoretical  mean  errors  are  indicated  by  solid  lines,  and  the  theoretical 
scatter  is  represented  by  twice  the  standard  deviation  of  the  error.  Figure 
10  shows  a comparison  between  the  predicted  mean  expected  coherence  and  the 
experimentally  derived  data. 

The  agreement  between  theory  and  experiment  is  remarkably  good  in  all  cases. 
Any  slight  errors  are  probably  due  to  inaccuracies  in  the  damping  factors  which 
were  derived  from  optimum  measurements  of  the  JdB  resonance  bandwidth. 


4.  APPLICATIONS  TO  MEASUREMENT  AND  PREDICTIVE  PROCEDURES  FOR  SDOF  SYSTEMS. 

In  order  to  apply  the  data  shown  in  Figures  1-5  it  is  necessary  to  define 
the  errors  which  are  acceptable  in  a particular  analysis.  Given  this  information, 
the  maximum  tolerable  analysis  bandwidth  may  be  read  directly  from  whichever  of 
Figures  1-5  is  appropriate.  In  addition  it  is  necessary  to  consider  three  other 
features  of  structural  analysis  techniques,  namely  the  effect  of  data  windows 
as  applied  in  digital  signal  processing,  the  application  of  frequency  translation 
to  measurements  of  lightly  damped  systems  and  the  non-uniform  accuracy  of  constant 
bandwidth  analysis. 
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4.1  The  effect  of  d&ta  windows  in  finite  bandwidth  analysis. 

In  general,  digital  signal  processing  techniques  require  that  a 
smoothing  window  should  be  applied  to  all  types  of  signals  with  the 
exception  of  signals  that  may  be  classified  as  being  "transient"  in 
nature.  The  effect  of  such  a smoothing  window  is  to  increase  the  nominal 
analysis  bandwidth  by  a factor  which  may  be  as  much  as  2.  Further  details 
of  the  precise  effects  of  particular  windows  may  be  found  in  (j),  but  in 
general,  the  use  of  any  data  window  produces  greater  inaccuracies  in 
measured  frequency  response  data  as  a result  of  the  implied  increase  in 
effective  analysis  bandwidth. 

4.2  The  application  of  frequency  translation  to  measurements  of 
ii Truly  damped  systems. 

If  an  N point  Discrete  Fourier  Transform  is  used  to  process  time 
series  data  in  the  frequency  range  0.<£<f  , then  the  hominal  bandwidth 
Af  is  given  by 

Af  = 2f 

o 

N 

The  application  of  a data  window  may,  however,  lead  to  the  effective 

bandwidth  Af  being  as  much  as  2Af.  Thus 
a 

Af  = 4f 

a o 

N 

If  the  Forlifer  criterion  is  applied  to  measurements  of  a lightly 
damped  system  of  resonant  frequency  fp  and  damping  factor  ^ , it  is 
required  that 


Afa  < 0.25 

For  ^ = 0.001,  Af  4 0.0005  fD 

an 


If,  for  example,  fD  = 0.5f  , the  Forlifer  condition  may  only  be 

satisfied  if  a 20,000  point  transform  is  used.  If  such  a facility  is  not 
available  the  only  solution  is  to  translate  frequencies  in  the  neighbour- 
hood of  fD  down  to  a baseband  defined  by  O^f^f  . In  this  way,  the 

required  resolution  may  be  achieved,  but  only  in  the  much  narrower  band 
defined  by 


fR  “ fa^f«fR 


(fa  + fb 


= V 


If  N is  subject  to  an  upper  limit  of  1024,  frequency  translation  is 
likely  to  be  required  whenever  the  damping  factor  is  less  than  or  equal  to 
0.02,  in  order  that  the  Forlifer  criterion  should  be  satisfied. 
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4 . 3 The  non-uniform  accuracy  of  constant  bandwidth  analysis. 

Since  the  3dB  bandwidth  of  a SDOF  system  resonance  increases 
linearly  with  the  resonant  frequency,  constant  accuracy  can  be  maintained 
over  the  whole  frequency  range  of  a multiple  degree  of  freedom  system, 
for  which  the  damping  factor  is  independent ‘.pf  frequency,  if  the  analysis 
bandwidth  also  increases  linearly  with  frequency.  This  may  be  achieved 
by  utilising  constant  percentage  bandwidth  analysis  in  which  the  centre 
frequencies  of  adjacent  bands  increase  logarithmically  with  frequency. 

For  discrete  frequency  analysis,  consistent  accuracy  may  be  similarly 
achieved  by  using  logarithmically  spaced  excitation  frequencies.  The 
variable  accuracy  of  constant  bandwidth  analysis  can  thus  be  an  important 
drawback  to  the  use  of  FFT  processing  since  accuracy  at  the  low  frequencies 
is  sacrificed  with  the  effect  that  the  high  frequency  portion  of  the 
spectrum  can  become  needlessly  over-resolved. 

5.  CONCLUSIONS 

If  Af  is  defined  as  the  analysis  bandwidth,  and  Af  as  the  3dB  bandwidth 

3.  S 

of  a system  resonance,  then  from  the  preceding  text  it  is  apparent  that  the 
application  of  the  Forlifer  condition  given  by 

Af  ^ 0.25 
&f 

s 

produces  minimal  errors  of  order  O.ldB  in  power  spectral  estimates  of  single 
degree  of  freedom  system  behaviour.  However,  it  is  also  apparent  that 
simultaneous  measurements  of  transfer  function  data  are  subject  to  much  greater 
errors  of  the  order  0.3dB  in  amplitude  and  as  much  as  14°  in  phase.  These 
phase  errors  may  be  particularly  important  if  further  data  reduction  by  modal 
analysis  is  involved.  In  these  circumstances,  it  may  be  necessary  to  tighten 
*he  Forlifer  condition  to  derive  six  estimates  within  the  3dB  structural 
bandwidth.  Thus 


Af  ^ 0.1? 

a ^ 

TT 

B 

In  addition  the  use  of  data  smoothing  windows  and  frequency  translation 
have  been  identified  as  important  contributory  factors  in  the  determination 
of  the  accuracy  of  measurements  and  predictions  of  transfer  function  data. 

Finally,  constant  percentage  bandwidth  analysis  may  be  recommended  since 
this  technique  retains  constant  accuracy  throughout  the  whole  frequency  range. 


A W Walker  (HSO) 

AWW/APG 
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FIG.  6 DEFINITION  OF  ELECTRICAL  SINGLE  DEGREE  OF 
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: 

J 


. 
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COMPARISON  OF  THEORETICAL  MEAN  ERRORS  IN  POWER  SPECTRUM  RATIC 
EXPERIMENTAL  MEASUREMENTS  OF  BEHAVIOUR  OF  TWO  SDOF  SYSTEMS 
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